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SINGULARITY CONFINEMENT AND FULL-DEAUTONOMISATION: 
A DISCRETE INTEGRABILITY CRITERION 

B. GRAMMATICOSi, A. RAMANF, R. WILLOX^, T. MASE^ AND J. SATSUMA^ 


Abstract. We present a new approach to singularity confinement which makes it an efficient 
and reliable discrete integrability detector. Our method is based on the full-deautonomisation 
procedure, which consists in analysing non-autonomous extensions of a given discrete system ob¬ 
tained by adding terms that are initially absent, but whose presence does not alter the singularity 
pattern. A justification for this approach is given through an algebro-geometric analysis. We also 
introduce the notions of early and late confinement. While the former is a confinement that may 
exist already for the autonomous system, the latter corresponds to a singularity pattern longer 
than that of the autonomous case. Late confinement will be shown to play an important role 
in the singularity analysis of systems with non-trivial gauge freedom, for which the existence of 
an undetected gauge in conjunction with a sketchy analysis, might lead to erroneous conclusions 
as to their integrability. An algebro-geometric analysis of the role of late confinement in this 
context is also offered. This novel type of singularity confinement analysis will be shown to allow 
for the exact calculation of the algebraic entropy of a given mapping. 


PACS numbers: 02.30.Ik, 05.45.Yv 


1. Introduction 

Singularity confinement [T] is a discrete analogue of the Painleve property [2] of ordinary dif¬ 
ferential equations, which infers the integrability of a given equation from the local structure of 
its singularities. The crucial requirement there is that singularities, the position of which depends 
upon the initial conditions, do not introduce multivaluedness (which in general makes it impossible 
to represent the solution of the differential equation as a function). Analogously, the singularity 
confinement approach is based on the local study of the singularities that appear in a discrete 
system. Here as well we are interested in singularities with positions that depend on the initial 
conditions of the system and singularity confinement requires those singularities to disappear after 
a few iteration steps, lest they lead to indeterminacies that make the construction of the solution 
of the system impossible. The relevance of singularity confinement as an integrability detector is 
strengthened by the fact that all discrete systems integrable through spectral methods, studied to 
date, have been shown to possess confined singularities. On the other hand, linearisable discrete 
systems in general do not satisfy the singularity confinement criterion [3] , in close parallel to what 
happens in the continuous case, where linearisable differential systems in general do not possess 
the Painleve property either [J]. 

If this parallel between singularity confinement and the Painleve property had been perfect, 
we would of course have been in possession of an efficient and convenient discrete integrability 
detector. However, the discovery of non-integrable systems with confined singularities called the 
usefulness of singularity confinement as an integrability criterion into question. The best-known 
example of such a mapping is the one proposed by Hietarinta and Viallet in , which we refer to 
as the H-V mapping: 

^n+l “1“ ^n—1 — “b (I) 

The pattern of its singularities is {xn+i = 0,x„+2 = oo,x„+3 = oo,x„+4 = 0} and, since Xn +5 = 
Xn , its singularity is confined. However the authors of [5] have shown numerically that this mapping 
exhibits large scale chaos and thus cannot be expected to be integrable. This of course raises the 
question what singularity confinement might mean in this case? Clearly, the confinement property 
is related to some subtle cancellations occurring when one iterates a rational mapping. These 
cancellations will in fact reduce the growth of the degree of the successive iterates. (We should 
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point out here that, as shown by Bellon and Viallet [6], while the degree itself is not invariant 
under coordinate changes its growth is invariant and thus characteristic for the mapping). In 
fact, when a mapping is integrable, these confinement-related cancellations slow down the degree 
growth to such an extent that, asymptotically, it becomes polynomial [7]. However, in the case 
of the H-V mapping, whereas some cancellations do take place, these do not suffice to curb the 
asymptotic degree growth which remains exponential. Such rapid growth is the signature of non- 
integrability. Hietarinta and Viallet [5] therefore introduced a quantitative measure of the degree 
growth of a rational mapping: its algebraic entropy. If d„ represents the homogeneous degree of 
the numerator or denominator of Xm the algebraic entropy of the mapping is given by the limit 
E = lim„_,.oo ^ log dn ■ For an integrable mapping the algebraic entropy must vanish. On the 
other hand, a non-zero value for £ implies exponential growth and is therefore an indication of 
non-integrability. In the case of mapping (P the algebraic entropy can be computed exactly [8] 
and is found to be £ = log( ^~^^'^ ). 

Curiously, this mapping remained essentially a singleton as far as counterexamples to the con¬ 
finement criterion were concerned (but see [5]), despite the fact that the authors of [S] presented 
general arguments for the existence of whole families of confining, non-integrable, mappings. The 
status of singularity confinement became even more complicated with the issue of late confinement 
[I0| . While standard practice in the implementation of singularity confinement, for mappings with 
parametric freedom, had been to enforce confinement at the very first possibility, it was not at all 
clear at the time why one should abide by this rule and why, for example, one could not postpone 
confinement until a later occasion. It turns out mm that when one opts for a late confinement, 
the resulting system will be non-integrable despite its singularities being confined (and despite the 
fact that when confinement is implemented normally, the resulting system might be integrable). 

These problems led to a certain distrust of singularity confinement as a method for detecting or 
deriving discrete integrable systems. Still, there has always existed a domain - for which we coined 
the term deautonomisation m - where this criterion continued to thrive and in fact furnished a 
slew of novel results. What we mean by deautonomisation, is to consider the free parameters of 
a mapping (which a priori take constant values) to be functions of the independent variable, the 
precise form of which has to be obtained through the use of a certain discrete integrability criterion. 
The rationale behind this approach lies in the relation the growth properties of the solution of a 
mapping bear to its integrability. In the deautonomisation procedure one starts from an integrable 
autonomous system, obeying the low-growth requirement, and one seeks to extend it to a non- 
autonomous form while keeping the same growth. In most practical applications however, the 
integrability criterion one uses is in fact singularity confinement. The reason being that, compared 
to techniques that rely on the calculation of the algebraic entropy, the confinement criterion has the 
immense advantage that one can examine each singularity separately, establishing the constraints 
on the parameters one at a time and not all at once in a hopelessly entangled way. 

It is precisely this very same deautonomisation approach that will be shown to reinstate singu¬ 
larity confinement as a reliable discrete integrability criterion. In |13j we introduced the so-called 
full-deautonomisation approach, and we claimed that this is the proper way to perform the sin¬ 
gularity analysis of a given mapping. If the system is integrable, the characteristic equations for 
the constraints that one obtains for the parameters, will only have roots with modulus I, whereas 
the presence of a root with modulus greater than I implies non-integrability. In what follows we 
shall first illustrate the power of this approach through several examples, after which shall give 
a detailed discussion of the problems that arise due to gauge freedom in the mapping and of the 
solution the concept of late confinement offers to this conundrum. 

2. The full-deautonomisation procedure 

The deautonomisation procedure consists in assuming that the parameters that appear in a 
mapping are functions of the independent variable and in using some integrability criterion, like 
singularity confinement, to fix their precise form. Standard practice when applying this procedure 
is to require that the (confined) singularity pattern of the autonomous mapping and that of its 
non-autonomous extension be identical. (An analogous requirement can be formulated whenever 
the algebraic entropy criterion is used: one then requires that the degree growths are the same 
for the autonomous and non-autonomous mappings). The full-deautonomisation procedure is an 
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extension of the standard one, where one introduces terms (with non-autonomous coefficients) that 
are absent in the initial mapping but which, when present, do not modify the singularity pattern. 
To put it in a naive but evocative way: one not only replaces Is by functions, but also selected Os. 

We shall illustrate this procedure and its implications for the integrability of a given mapping 
on two examples. The first one is a mapping of the form 

1 ^ ^ 
^n+l ^n—1 — (2) 


Its singularity pattern is {0,oo^,0}, where by oo^ we mean that if we introduce a small quantity 
e and assume that Xn is finite and Xn+i = £, then Xn +2 will be of order 1/e^. Deautonomising ([2]) 
then consists in replacing the numerator of the right-hand side by a function a„ and to require 
the mapping to have confined singularities with exactly the same pattern as the autonomous one. 
This yields the constraint Un+i = a-n-i, which gives an integrable, but trivial, non-autonomous 
extension of In fact, by introducing the appropriate gauge a;„ —>■ jnXn, with 7 )^ = a^/a„-i 
and 7n-i-i = 7n-i, we can put a„ = 1 for all n. 

In order to proceed to the full-deautonomisation of (I2|) we must add terms that do not modify 
the initial singularity pattern. It is straightforward to convince oneself that the only possible such 
extension is by adding a term inversely proportional to x, which leads to 


^n+l T ^n—1 



Xn 


(3) 


Requiring the initial singularity pattern to be confined ~ i.e. if we start from Xn+i = e, that Xn+i 
takes a finite value which depends on Xn - we find the constraints a„+i = a„_i and bn+i — 2 &„ - 1 - 
bn-i = 0. The integration of the latter leads to the characteristic equation (A — 1)^ = 0 and finally 
to bn = an + /3. This expression, combined with a„ = 1, turns (|31) into a well-known form of the 
discrete Painleve I equation [T^ . 

The second example is a non-integrable mapping with confined singularities of a form similar 
to ([ 2 ]), 

Xn+l T Xn—1 — T 7 (^) 

Xn 

but with singularity pattern {0,oo'^,0}. Its standard deautonomisation follows exactly that of ([U 
and leads to the same constraint, and the same uninteresting result. Its full-deautonomisation is 
however of particular interest. Including in (|4]) all the terms that leave the singularity pattern 
unchanged, leads to the mapping 


Xn+l + Xn-1 — -1- J H-^ H-J. 

^r> ^r> 


(5) 


The confinement constraints obtained by requiring that Xn+i be finite, just as in the integrable 
case above, are: a„+i = a„_i, dn+i = -d«-i, c„+i = c„_i and 6 „+i - 46„ -|- 6 „_i = 0. The 
first three equations can be integrated in terms of periodic functions, but they can also simply be 
satisfied by the choice a„ = 1 , = 0 , c„ = 1 or even c„ = 0 , without influencing the reasoning 

on the fourth constraint. Solving the latter, we find the characteristic equation A^ — 4A -I- 1 = 0 
for which, contrary to the integrable case above, the roots are not roots of unity but are equal 
to 2 ± v^- As we posited in [T3], the existence of a root in the characteristic equation for the 
coefficients, with modulus greater than 1, is an indication of non-integrability. Our argument is 
further strengthened by the observation that the logarithm of the largest root is precisely the value 
of the algebraic entropy that can be computed using the standard numerical procedure. 

Thus our claim is that the proper way to perform the singularity analysis using the singularity 
confinement criterion is by performing a full-deautonomisation of the mapping. While our argu¬ 
ments may appear heuristic at this stage it is possible to make them more rigorous by a detailed 
algebro-geometric analysis of these mappings. 


3. An ALGEBRO-GEOMETRIC JUSTIFIGATION OF FULL-DEAUTONOMISATION 

In this section we analyse the mappings ([31) and from the point of view of algebraic geom¬ 
etry. In particular, we shall derive conditions on the parameters, by blowing up the mappings at 
their indeterminate points, such that, ultimately, each mapping induces an automorphism on an 
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appropriately structured rational surface, obtained by blowing-up x P^. This requirement may 
well be stronger than the singularity confinement criterion itself (we plan to clarify this point in 
a future publication) but it will be shown that the conditions on the parameters we obtain from 
this analysis do coincide with those obtained in the previous section. Moreover, the role these 
conditions play in determining the algebraic entropy of a mapping will also become clear. 


3.1. The case of an integrable mapping. Let us start with ([3]), which we shall interpret as the 
following birational map on P^ x P^: 

tfn- P^XP^—^P^XP\ {Xn,yn) ^ iXn+l,yn+l) = (yn,-Xn + — + (6) 

V Vn yiJ 

We also introduce the variables Sn = l/x„ and = 1/yn and we cover P^ x P^ with four copies 
of C^, as: 

pi X P^ = {Xn, yn) U {Xn, tn) U (s„, ?/„) U (s„, t„). (7) 

Obviously, the point {sn,yn) = (0,0) is an indeterminate point for the map (pn which we shall try 
to regularise by a ‘blow-up’, i.e. by introducing two new coordinate charts that will replace the 
single point (s„, 2 /„) = (0, 0) (which will be called the base-point for the blow-up): 


{Sn,yn) ^ 


The map can now be expressed as 

1 

yn+l — 2 Sn 



( 8 ) 

(9) 


which, however, is now indeterminate at y^^ynj = (0,0) in the second new coordinate chart. 
Hence we introduce two new coordinate charts 





\yn’ Sn) 



( 10 ) 


at the base-point i^tyn) = (0,0), after which the mapping is given by 


yn+l — 2 s„ [ 1 + ^"2/" 2 + 

in terms of the second coordinate chart in However, this does not yet resolve the indetermi¬ 
nacy, which now manifests itself at ^p-, 0^, which also requires blowing-up. Note that 

this base-point now depends on the independent variable n through the parameter function a„, 
whereas the base-points in the previous two blow-ups were independent of n. This third blow-up 
is carried out by means of the coordinate charts 


Sn 

vl 






in terms of which we obtain 


yn+l — 




( 12 ) 


(13) 


which still has an indeterminacy at the point ^p — ^ j , j = Oj . Note that this point 

now also involves the second parameter, bn- A blow-up at this base-point, using the coordinate 
charts 



( 14 ) 
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finally lifts the indeterminacy, as one obtains: 


1 


Un+l — Sn 






(15) 


This sequence of four succesive blow-ups at the point {sn,yn) = (0,0) is depicted on the left 
in Figure [TJ together with the resulting configuration of exceptional curves. These curves can be 
thought of as (local) copies of P^, spliced in as it were, at the base-point of each blow-up and that 
are described by the coordinate charts introduced for that blow-up. 



tn = 


0 


■ yn = o 


Xn = 0 


I = 0 


Figure 1. Diagram showing, on the left, the base-points in the successive blow¬ 
ups of the map ipn and the exceptional curves resulting from these blow-ups. Same 
on the right for the indeterminate points of the inverse map 
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However, this does not regularise all indeterminate points for the mapping ipn- The easiest way 
to resolve the remaining indeterminacies and lift ipn to a well-defined automorphism, is to consider 
its inverse map 

ip~^: P^xP^--4 P^xP\ (a;„+i,y„+i) (x„,?/„) = (-y„+i-f-f) , (16) 

V Xn+l <+i J 

which is undefined at (xn+i,tn+i) = (0,0). As was the case for this indeterminacy can be 
resolved by four succesive blow-ups, the base-points and coordinate charts of which are given below 
(the resulting mappings are however omitted, for brevity). 

blow-up at (x„+i,t„+i) = (0,0) : (a:„+i, t„+i)-«- ( x„+i,) U ( ) (17) 

V Xn+lJ \tn+l J 


blow-up at 




^n+li I — (0; 0) ■ 




^n-|-l 5 9 

^n +1 


u 




(18) 


blow-up at 


blow-up at 





These blow-ups and their corresponding exceptional curves are shown on the right in Figure[T] Note 
that, here as well, the base-points in the last two blow-ups depend on n, through the parameters 
a„ and 6„. Moreover, it is easily verified that the full set of eight blow-ups described above resolves 
all indeterminacies in the autonomous case of (ED, i.e. when the parameters a„ and are constant 
for all n. In the non-autonomous case however, these parameters have to satisfy certain constraints 
for the blow-ups to fully resolve the indeterminacies in the mapping. Indeed, it is easy to check 
that ipn maps the point {xn,tn) = (0,0), i.e. the first indeterminate point for into the first 

indeterminate point {sn+i,yn+i) = (0,0) for ipn+i- Similarly, (pn maps the second indeterminate 
point (^Xn, j = (0,0) for into the second indeterminate point , yn+i^ = (0, 0) of Pn+i 

(in its (ED avatar). 

In order to analyse the relations between the other indeterminate points, we define: 


Pn ■■ 




Qn+1 



^ yn+1 


Pn ■■ 



, yn+1 




( 21 ) 

( 22 ) 


Rti 


I ft 


Xr, 


Xr 


1 


Xt, 


5rt+l 


^n—1 

1 / Srt+l _ 1 


= 0 ,— 


^n—1 


yn+1 


R„ 


\yl+i 

f ^n-hl 


<^n+l 


1 yn-\-l I — 




,0 




1 


yn+1 \yn+l 0,n+l ) 


!yn+1 I — 


“t” ^n—1 


(23) 

(24) 

(25) 


and 
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It is easily verified that Pn is the image of Pn under the mapping cpn- Moreover, in the au¬ 
tonomous case (a„ = a, bn = b,'^n) it is clear that Pn = Qn+i, as depicted in Figure [5] 



Figure 2. The link between the base-point Pn for the third blow-up for ^Pn-i 
and the base-point Qn+i for the third blow-up for i^n+i- 

The same is true for the points = ^Pn{Rn) and .Sn-i-i which automatically coincide in the 
autonomous case (cf. Figure |21)- Hence, in the autonomous case, since (pn maps the indeterminate 
points of Pn-i to those of Pn+i (and vice versa, mutatis mutandis) and since all these indetermi- 
nacies have been resolved, in general, by blowing up x P^ eight times at every value of n , the 
mapping pn has been successfully regularised. Now, it is clear that these same eight blow-ups of 
P^ X P^ can also fully regularise the non-autonomous mapping if one imposes certain constraints 
on the parameter functions. In particular, and Qn+i will coincide if and only if a„+i = a„_i 
holds, which is nothing but the first constraint obtained from singularity confinement for the map¬ 
ping ©• Furthermore, under this constraint we can also require = Sn+i, which amounts to 
requiring bn to satisfy bn+i — ^bn + bn-i = 0, exactly as in the singularity confinement analysis. 



Figure 3. The link between the base-point Rn for the fourth blow-up for Pn-i 
and the base-point 5'n-i-i for the fourth blow-up for Pn+i- 

It is also clear that taking bn = 0 (''n) in no way affects the number of blow-ups needed to fully 
regularise the mapping, nor does it influence the configuration of exceptional curves obtained from 
these blow-ups. This statement is in fact equivalent to saying that the presence or absence of the 
term bn/xn in the right-hand side of the mapping ([3]) does not alter its singularity pattern, which 
is the fundamental observation that underlies the full-deautonomisation procedure. To explain the 
precise link between the singularity pattern and the blow-ups carried out above, it is important 
however to take a closer look at the configuration of exceptional curves that is obtained from them. 

The sequence of blow-ups of x P^ described above, produces a rational surface [13] that 
contains all the exceptional curves obtained from the blow-ups (these consist of all the points that 
correspond to the base-points of the blow-ups) and which is obtained by glueing together all the 
different coordinate charts used in the blow-ups. In fact, to be completely accurate, one should 
say a family of surfaces, since the base-points for four of the blow-ups are n-dependent and the 
corresponding surfaces will therefore not be strictly identical. However, algebraically speaking, 
only the relative positions and, especially, the intersection pattern of these curves matters and one 
can therefore represent this entire family of surfaces by a single diagram, such as in Figure |31 The 
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curves depicted in this figure are part of a free Abelian group, the so called Picard group [M], which 
is in this case of rank 10 as it is obtained by 8 blow-ups from x P^. In fact, we can take the curves 
£> 1 , £> 2 , ■ ■ ■, Ds,Ci and C 2 , as a basis that will generate the entire group. Note that the curves 
labelled with D differ in nature from the curves Ci and C2 as the latter have self-intersection -1 
whereas the former have self-intersection -2 and, hence, can not be trivially ‘blown-down’. In fact, 
the intersection pattern of the curves with self-intersection -2 has the form of a Dynkin diagram of 
type £Iy^^. This type of characterization of the surface plays an important role in the classification 
of discrete Painleve equations obtained by Sakai where the non-autonomous mapping ([3]) that 

satisfies the two conditions a „+2 = a™ and 6„+2 — -|- &„ = 0 is classified as a discrete Painleve 

with symmetries associated with the affine Weyl-type group 

Di 



Figure 4. Diagrammatic representation of the family of surfaces on which ipn 
acts as an automorphism. The labels {x = 0} and {y = 0} refer to the strict 
transforms of the corresponding curves in P^ x P^. 


Moreover, as is clear from the correspondence between the indeterminate points of ip~^ and 
ifn+i under the mapping explained above, ipn induces the following automorphism on the 
Picard group: 

Di —>■ £>2 —^ D 3 —^ £>4 —>■ Di , £>5 —^ £>6 —^ D 5 , Dy —>■ £)§ —^ Dy (26) 

{y = 0} —> Cl —>■ C 2 —> {x = 0} , (27) 

where {a: = 0} is found to be linearly equivalent to 

-Di+ 2£»2 - £>3 + £>4 - -D 5 + 2£>6 + Ds-Ci+ 2 C 2 . (28) 

Note that the sequence of values for y in (EH) is nothing but the singularity pattern {0, 00 ^, 0}. 

With respect to the basis for the Picard group we have chosen, this automorphism, which we 
shall denote by can be written in matrix form as: 


0 0 0 1 

0 -1 

10 0 0 

0 2 

0 10 0 

0 -1 

0 0 10 

0 1 

0 1 

0 -1 

1 0 

0 2 

0 1 

0 0 

1 0 

0 1 

0 

A 


(29) 


v 


/ 
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where . It is well-known [5] that the algebraic entropy of a birational map such 

as can be obtained from the largest eigenvalue, A*, for the induced automorphism 

8 ■.= lim — log(deg(v?„+m-i‘/7n+m-2---</5„)) =log|A*|, (30) 

m—>-oo 777, 

where by deg((^), we mean the maximum of the degrees of the numerator and denominator of 
the birational map ip. Of course, as is clear from the above matrix, its upper-left sub-matrix 
(corresponding to the sequences (l26ll of D's) cannot yield an eigenvalue with modulus greater than 
1 since it is unitary. Hence, the part of the matrix that decides on the positivity of the entropy, 
i.e. on the integrability or non-integrability of the mapping (/?„, is the sub-matrix A which has 

det(A/-H) = A^ - 2A-f 1, (31) 

as characteristic polynomial. All eigenvalues of therefore have modulus 1 and the algebraic 
entropy of tpn is zero. The mapping, subject to the conditions a „+2 = On and bn +2 — ‘2bn+i + bn = 0 
which are needed in the above regularisation, is therefore integrable. It is important to note that the 
characteristic polynomial for the A part of ip^ is identical to that for the only non-trivial constraint 
on the parameters in the mapping bn +2 — ‘2-bn+i + bn = 0 (as an such that an +2 = a-n, 
corresponds to a simple gauge transformation and, as such, has to be considered as a trivial 
constraint vis-a-vis the integrability of the mapping). We shall see on a second example that 
similar conclusions can be drawn in the case of non-integrable mappings. 

3.2. The case of a non-integrable mapping. Next, we consider the special case of mapping 
when Cn = dn — 0. The analysis that follows can of course also be carried out for the general case, 
i.e. for non-zero c„ and dn, but besides complicating the relevant mathematical expressions to the 
point where they become almost illegible, such a generalization does not alter the conclusions we 
shall obtain in the restricted case. 

Here as well, we interpret the mapping as a birational map on x P^: 

P^XP^—^P^XP\ iXn,yn) ^ iXn+l,yn+l) = (yn,-Xn + — + (32) 

V Vn ytJ 

which is undefined at {sn,yn) = (0,0) (in the coordinates introduced in ([7])) and eight successive 
blow-ups are needed to lift this indeterminacy. However, as their precise mathematical expres¬ 
sions become quite involved, we prefer to give these 8 base-points (i = I-- - ,8), with the 
corresponding coordinate charts and the maps obtained after each blow-up, in an appendix to this 
paper. 

As in the previous case, the remaining indeterminacies are best resolved at the level of the 
inverse map 

P^xP^-->P^xP\ (a:„+i,y„+i) !-)■ (a;„,?/„) = (-y„+i-f-f] , (33) 

V Xn+l 4+1 / 

which is indeterminate at the point (a:„+i, t„+i) = (0, 0). Again, eight blow-ups are required to lift 
the indeterminacy and the base-points Pn\i (* = 1 • • • , 8) and coordinate charts for these blow-ups 
are also given in the appendix. 

Now, it can be easily checked that </)„, maps the base-points of the first four blow-ups for 
ipn-i into the base-points Qn+i of the corresponding blow-ups for (j)n+i, i-e.: 4>n = Qn+i 

for i = 1, 2,3 and 4. For the base-points of the fifth blow-up we find that 

matches up automatically with the base-point Qn+i ■ , 2/n-i-i^ = for the fifth blow¬ 

up of (jjn+i, in the autonomous case, but that the parameter a„ needs to satisfy the constraint 

Uti+I — Un—l, (i^^) 

for this match to occur in general, in the non-autonomous case. Under this constraint, the images 
under (/)„ of the next two base-points for 4>n-i again coincide naturally with those for 4>n+i, 
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i.e.: (j}n = Qn+l^ ior i — 6 and 7. However, the remaining base-points only match up 

automatically in the autonomous case. For the general non-autonomous mapping, one finds that 


will coincide with 


1 


1 


^,3 

yn+1 


/ ^n+l 
XVn+l <^n+l 


^ iJn+l I — 


— 4bn + bn-1 


^n +1 ■ 


1 


/ 5w+l 


1 


Vn+l \2^ra+l fln+l 


) Vn+l I — 


^n -1 


'Jn+1 


^n +1 


if and only if, under the constraint (jd5|) . satisfies the condition 

^n+i 4bn -t- bji—i = 0. 


,0 


(36) 


(37) 


(38) 


We have thus recovered the two conditions obtained from the singularity confinement analysis 
of the mapping ([S]) in section [21 and we can conclude that if the parameters a„ and 6 „ satisfy the 
constraints the map (j)n can be fully regularised through (two times eight) 16 blow-ups. 

In other words, under the above conditions, (/>„ will act as an automorphism on the (family of) 
surfaces obtained by glueing together all the coordinate charts used in the blow-ups. Figure 
depicts the exceptional curves obtained from the blow-ups, that are contained in such a surface. 

Note that, as in the first example, here as well we see that the presence or absence of the term 
bnixn in the mapping does not alter the number of blow-ups (or their structure) needed for the 
regularisation of the mapping, which was the basic ansatz in the full-deautonomisation procedure. 


Di 



Figure 5. Representation of the exceptional curves in the surface on which 
(EH) acts as an automorphism. The labels {x = 0} and {y = 0} refer to the strict 
transforms of the corresponding curves in x P^. 

Since 16 blow-ups of P^ x P^ were needed to construct the surface shown in figure 0 its Picard 
group will be of rank 18 and we can take the curves Hi, II 2 , ■ • ■, Hig, Ci and C 2 as a basis. Moreover, 
the analysis of the correspondences between the base-points for the different blow-ups, given above, 
also establishes the automorphism induced by the map </)„ on the Picard group: 

Di —^ D 2 —>■ D 3 —>■ D 4 —>■ Di , H 5 —>■ Dq —>■ Dj —>■ Hg —>■ H 5 , 


Hg —^ Hio —^ Hii —^ Hi 2 —^ Hg , 


Hi3 —Hi4 —H 13 , 


Hi5 Hi6 —Hi5 , 


(39) 

(40) 
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and 

{ 1 / = 0} —>■ Cl —>■ (72 —^ {a; = 0}. (41) 

Note that the y dependence in this last sequence of curves corresponds exactly to the singularity 
pattern {0, 00 ^, 0} for the mapping ® . 

As {cc = 0} is linearly equivalent to 


— Di + 4 i 72 — 7?3 + -D 4 — Dq + 4:Dq — ZI 7 + 2i7g — Dg + 4i7io — Du + 3Z?i2 

~ Di3 + 4i7i4 + Dig — Cl + 4 C 2 , (42) 


( O' * \ 

— 0 —I, where ct is a unitary 


0 -1 


sub-matrix (as it represents the permutations (I39I40I) of the D curves) and where A = ^ 

Hence, it is the characteristic polynomial for A, 

det(A/-A) = - 4A-b 1, (43) 


that determines the largest eigenvalue of : A, =2 + This proves that the algebraic entropy of 
(pn - and therefore also that of ([S]) - is positive and that this mapping is indeed non-integrable, as 
was surmised in section [2J The important thing to note however is that, here again, the character¬ 
istic polynomial for the sub-matrix A coincides with that of the (only) non-trivial constraint, (13811 . 
on the parameters in the mapping. Because this is exactly the factor in the full characteristic poly¬ 
nomial of (j)* that decides on the non-integrability of the mapping, it becomes clear that although 
the absence or presence of the bn/xn term in the mapping does not alter the surface obtained from 
the blow-ups - or the action of the mapping on the Picard group for that matter - without this 
term, this correspondence between the constraints on the parameters and the action on the Picard 
group would be broken and it would be impossible to read-off the value of the largest eigenvalue 
of (j)* just from the constraints on the coefficients. This fact justifies the full-deautonomisation 
approach, in which one supplements the mapping with terms that do not alter the singularity 
pattern, but that are general enough to capture, so to speak, the full dynamics of the mapping 
induced on the Picard group, at the level of the parameters. 


4. Applications of the full-deautonomisation procedure 

In the preceding sections we presented the basic workings of the full-deautonomisation method 
and its justification based on algebro-geometric arguments. Here, in order to underline the new 
power the singularity confinement method gains from this approach, we present selected examples 
where this novel method makes it possible to detect non-integrability despite the fact that they all 
have confined singularities. 

In [S], Viallet offered an example of a mapping that possesses confined singularities but which 
has non-zero algebraic entropy and is thus not expected to be integrable: 

1 , , 
Xn+lXn-l = Xn-\ -h U (44) 

Xn 

(The above mapping, with a = 0, was first published by Hone [16] who, however, attributes its 
paternity to Viallet). The singularity pattern of (|T]) is {0, oo, oo^, oo, 0}. The standard deau- 
tonomisation of (PI) consists in just assuming that a„ is a function of the independent variable. 
Implementing the singularity analysis in this case, while requiring that the pattern be the same 
as for the autonomous one (and in particular when starting from a finite and a vanishing Xn+i 
requiring that Xn +7 be hnite), we obtain the constraint Qn+e = On, he. a„ should be a periodic 
function with period 6. 

The full-deautonomisation of (1441) would require that we introduce functions multiplying the 
terms proportional to x and 1/x in the right-hand side as well as keeping a„. However it turns 
out that a does not play any role in fixing the value of the algebraic entropy, and since a„ = 0 is 
a possible solution to the constraint that was obtained in the previous paragraph, in what follows, 
we choose its value to be 0. We shall thus work with 

_ 1 I 

Xn+lXn—l — bjiXji -\- 


( 45 ) 
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We perform the singularity analysis requiring that the singularity pattern be the same as for the 
autonomous case and obtain the confinement constraint 

bn+ebnC^+^C^+l = ^„+4&„+3^n+2Cn+6Cn- (46) 

Inspecting equation (H51) it is clear that there exists a global gauge freedom which allows to put 
either or c„ equal to 1. Thus the confinement constraint can be expressed as either 

'^n+5^n+l ~ Cn+eCn (47) 


or 

bn+ebn = 6^+4fe^+3^n+2- (48) 

Seeking a solution to these constraints for log c„ or log in the form A" we obtain the characteristic 
equation (A^ ± A + 1)(A'^ — A^ — 2A^ — A + 1) = 0, where the plus and minus signs in the first 
factor correspond to gB and (gSll respectively. The largest root of this characteristic equation is 

(1 + •\/T7)/4 + Y^(l + •\/T7)/8 the approximate value of which is 2.081019, in perfect agreement 
with the value obtained by Viallet in his calculation of the degree growth of (l4^ . 

Next we turn to a generalisation of mapping (g]) which was analysed in section gj 

^n+l 4“ ^n—1 — fT ^ (49) 


where /c is a positive integer greater than 2 (the cases k = 1 and k = 2 being integrable). The 
requirement that the singularity confines after three steps, introduces the constraint a„+i = 
(—l)^a„_i. (For even k this means that Un must be a periodic function with period 2 while 
in the case of odd k the function a„ must be of period 4 but involving only powers of i and —i). 
In order to show the non-integrability of (I49L we perform a full-deautonomisation by adding all 
terms of the form x~^ with £ = 1, - ■ ■ , fc — 1 with free coefficients. However, with hindsight, as in 
the case analysed in section [21 the term x~^ alone already suffices to obtain the answer we seek. 
We thus work with 


^n+l 4 “ ^n—1 



X„ 


(50) 


and we require that the singularity pattern be { 0 , oo*, 0 }, as it was in the absence of the bn/xn term. 
Using the proper choice for we readily find the condition — kbn + bn-i = 0 which leads to 
the characteristic equation A^ — fcA 4 - 1 = 0, the largest root of which is equal to (fc 4- y/k'^ — 4)/2. 
A numerical estimate of the degree growth in the case fc = 4, obtained with the help of the 
Diophantine approximation introduced by Halburd leads after 12 iterations to a value of 
3.732052, very close to the expected one, 2 + v^. 

The following variant of (I49|) is also interesting in the sense that it corresponds to a different 
singularity pattern. Indeed for 

Xn-\-l 4“ Xji—i — 1 4“ (61) 

Xr) 


and fc > 1 , we find the pattern { 0 , oo^, 1 , oo^, 0 } provided a„ satisfies the constraint a„+ 24 -a „_2 = 0 . 
Here again we proceed to the full-deautonomisation by adding a term inversely proportional to x 
with a free coefficient: 


Xn+i 4- a:„_i — 1 4-1— jr- (52) 

Xn X^ 

After implementing the constraint for we find that, for the singularity to be confined, must 
satisfy: 6„+2 — k{bn+i + 5„_i) -|- bn -2 = 0, with characteristic equation A^ — fc(A^ 4- A) 4- 1 = 0. 

The largest root for the latter is (fc 4- Vk'^ 4 - 8)/4 -|- {k'/W+ 8 4 - — 4)/8, which is greater 

than 1 for A: > 1. Thus we expect m to be non-integrable for fc > 1. (The case A: = 1 is 
a well-known integrable one, obtained from (1521) with a„ = 0 and satisfying the constraint 

bn+2 — bn+l — bn-l 4- &n-2 = 0) [T]. 

The final example of this section is an extension of a mapping presented in m 


Xn+lXn-l = (53) 

xt + 1 

where k is an integer greater than 2. (The cases A: = 1 and k = 2 are known to be integrable) [T5] . 
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The cases of even and odd k should be treated separately since their singularity structure is 
different. In the case of even k the singularity patterns are very simple. We have {r, 0, —1/r}, 
where r is a fcth root of 1, and {s,c», 1/s} where s is a fcth root of —1. For odd k the singularity 
patterns are {r, 0, —1/r, oo, —r, 0, l/rj and {s, oo, 1/s} where, again, r and s are fcth roots of 1 and 
— 1, respectively. Just as in the case fc = 4 analysed in m. the full-deautonomisation of can 
be limited to the introduction of a single free function: 


_ _ _ *^71 

^n+l^n—1 — t , 1 
^n + 1 


(54) 


In the case of even k we analyse the singularity Xn = a™ and require that it confine in three 
steps, exactly as in the autonomous case. We obtain a confinement condition a„+ia„_i = 
which leads to the characteristic equation — fcA + 1 = 0. The largest root of the latter is 
{k + \/fc^ — 4)/2 which for /c > 2 is always greater than 1. The same singularity for odd k leads 
to a confinement constraint a„+ 3 a „_3 = a^+ 2 ®n- 2 - The corresponding characteristic equation is 
now A® — fc(A^ + A) + 1 = 0. While it does not appear possible to express the largest root of 
this equation in terms of radicals, it is easy to obtain the first terms of its asymptotic expansion, 
k + l/k^ — 1/k^, which offers a good precision already in the case fc = 3. Thus for all k greater than 
2 the singularity confinement approach, combined with the full-deautonomisation procedure, offers 
a clear indication of the non-integrability of (j53|) . despite the presence of confined singularities. 


5. Early versus late confinement 

In m we examined a particular application of singularity confinement, known as late confine¬ 
ment. This term however requires some clarification. When applying the singularity confinement 
criterion in conjunction with some deautonomisation approach, it is standard practice to confine at 
the very first confinement opportunity. The notion of “very first” opportunity, however, becomes 
hazy when the initial mapping has several singularity patterns, leading to more than one possible 
deautonomisation. This is something that can be dealt with, though it requires some experience, 
especially as none of these deautonomisations are what we call “late” ones, although some of them 
correspond to singularity patterns that are longer than others. What can serve here as a guide 
for the perplexed is that one has to require that every singularity pattern for the non-autonomous 
extension be identical to one of the autonomous mapping, but in the non-autonomous case this is 
not the only possibility. As was first pointed out by Hietarinta and Viallet [10], it is possible to 
bypass the first confinement opportunity and still impose confinement at some later stage. However 
the consequence of such a “late” confinement is that it leads to a non-integrable system while a 
normal confinement would have led to an integrable one. The algebro-geometric analysis of this 
phenomenon was presented in m- 

Another possibility is that a confinement possibility appears earlier in the singularity pattern. 
The existence of such a possibility is something that can be assessed already in the autonomous 
case. While the late confinement situation is always possible (something that is obvious once one 
follows the algebro-geometric analysis), an “early” confinement possibility is a particular situation 
which may arise or not, depending on the mapping. Still, several examples do exist and in what 
follows we shall illustrate the consequences of early and late confinement on an otherwise integrable 
mapping. 

We start with the mapping 

= a„(I - a;„). (55) 

In the autonomous case, i.e. when o„ = a for all n, we have the singularity pattern 

{1, 0, a, oo, oo, a, 0,1}, (56) 

meaning that if we start from a finite Xn and Xn+i = 1, we find that Xn+g is finite and in fact equal 
to Xn ■ The standard deautonomisation of (1551) was presented in [T5]. The singularity pattern, with 
the same structure as for the autonomous case, is now 


}1, 0, ang-2: OO, OO, Cn+gan+A/a.n-\-2 ^ Hi ajig-7a.n-\-2/ (an-t-5^n-t-4)}- 


(57) 
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The confinement constraint is a„+ 7 a „+2 = an+ 5 an +4 and can be integrated to log a„ = an + /3 + 
7 (— 1)" + Sj^ + where j = e~^ . Under this deautonomisation equation, (l55ll is a discrete 

Painleve equation associated with the affine Weyl group . 

However, there exists another possibility. Instead of confining after 8 steps we can choose to 
confine just after the third step, by taking a„ = 1. This confinement is indeed possible but 
its consequence on the mapping is drastic: the mapping is not only autonomous but, starting 
from Xn and x„+i we find the succession of values (1 — x„+i)/x„, (x„ + x„+i — l)/(x„x„+i), (1 — 
Xra)/xn+i, x„, x„+i, • • •, i.e. the mapping is now periodic with period 5. In fact, in all studied cases 
we have observed the same phenomenon: if early confinement is enforced, the mapping becomes 
periodic (unless it becomes entirely trivial). 

We turn now to the case of late confinement by deciding not to confine after 8 steps. In this 
case Xn+9 is infinite and the next confinement opportunity appears at the level of Xn+iA- The 
confinement constraint is now a„+ioan+ 5 a„ = a„+ 8 an+ 7 an+ 3 a„+ 2 ) the characteristic equation for 
which is — A® — + A® — + 1 = 0, with largest root A = 1.293486. A numerical 

evaluation of this root using Halburd’s Diophantine method yields, after 50 iterations, a value of 
1.31 and the convergence to the value obtained from the characteristic polynomial is extremely 
slow. In the next section we shall prove however that the above largest root does indeed lead to 
the correct degree growth for this mapping. 

One further remark should be made at this point. The late confinement we introduced at the 
level of x„+i 4 is not the only possible one. In fact, infinitely many confinement opportunities 
exist every 5 steps, i.e. for all indices of the form n + 9 + 5m. A study of the algebraic entropy 
for successive late confinements shows that it grows monotonically and approaches an upper limit, 
which can be computed exactly. Starting from the characteristic polynomial for the first, standard, 
confinement, Po(A) = A® — A® — A^ + 1, it is easy to express the successive Pm in terms of Pq as 

Pm{X) = (Po(A) - l) + 1. (58) 

Since the largest root of Pm has modulus greater than I, we find, at the limit m —>■ oo, that 
this root must coincide with that of Po{X) — 1, which is given by the equation A® — A — I = 0. 

We obtain thus Aoo = ^(i/27/4 — y^23/4)®/® + (•\/27/4 — y^23/4)“^/®^ /v^ (with a numerical 
value approximately equal to 1.324718), the logarithm of which constitutes the upper limit of the 
algebraic entropy of the system, one that would necessitate an infinitely delayed confinement. 

6. On the pitfalls of gauge freedom 

The analysis in the previous section was based on the special form (IHKll of the mapping, 
which corresponds to a precise choice of gauge. However, another gauge could for example give 
x„+ix„_i = 1 — c„x„, and in general, for a given mapping, it is not always clear what the optimal 
choice for the gauge might be. Moreover, in what follows we shall explain that a ‘bad’ choice of 
gauge, combined with an inadvertent analysis of the constraints obtained from singularity confine¬ 
ment, may actually lead to entirely wrong conclusions. As we shall see however, this problem is 
by no means insurmountable. 

In order to illustrate this subtle point, we start from the mapping (1551) and use the full gauge 
freedom, which allows us to write it as 

X„+ 1 X „_1 = an- bnXn- (59) 

Next we repeat the singularity analysis starting from a finite Xn and Xn+i = On+i/^ri+i, while 
requiring confinement after 8 steps. We readily obtain the condition 

an-|-8Qn+l^n+5&ra+4 = On+7an-|-2feri+8^n+l ■ (60) 

The case analysed in the previous section corresponds to = a„, but any relation between a„ 
and bn is a priori conceivable. We can take for instance bn = a^, for some power k ^ 1, in which 
case the characteristic polynomial obtained from the confinement condition becomes (fc — 1)(A^ + 
1) -I- A® -I- A — fc(A^ -I- A®), which factorises into (A — 1)^(A -I- 1)(A^ -I- A -I- 1) [(/c — 1)(A^ -I- 1) -I- A]. 
While the first three factors correspond to the integrable n-dependence obtained in the previous 
section, the last factor creates a major problem. In fact, choosing k appropriately, we can make 
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the largest root of this characteristic polynomial equal to any positive number we wish. In which 
case, according to the conjecture presented in |13] and which is supported by all the results in the 
preceding sections, this would be an indication of non-integrability for (1591) . This conclusion, if 
correct, would of course toll the bell for singularity confinement as a discrete integrability detector, 
since it is in direct contradiction with the basic intuition that says that the mapping (ESI) should 
be integrable, regardless of the gauge one chooses. 

However, this hasty conclusion is totally unwarranted. Using the transformation Xn —t "fnXn, the 
mapping (15^ with can be cast into the form (ISSl) by choosing the gauge function 7 „ such 

that it satisfies 7 „ = ( 7 „+i 7 „_i)^“^, which of course corresponds to the characteristic equation 
{k — 1)(A^ + 1) + A = 0, for log 7 „. Hence it becomes clear that the troublesome factor obtained in 
the above analysis is entirely due to the specific choice of gauge in the mapping dMl), and that its 
contribution to the coefficients can be easily undone by a simple transformation of the dependent 
variables. Most importantly, this factor in the characteristic polynomial does not influence the 
degree growth of the mapping at all, as we shall explain shortly by means of an algebro-geometric 
analysis of dMI). Moreover we shall show that, strictly speaking, such an analysis is not even 
necessary to solve this gauge conundrum. 

In fact, as will become clear from the algebro-geometric analysis of dMl), the key to successfully 
detecting contributions to the confinement constraints that are due solely to gauge freedom, lies in 
the late confinements of the mapping. Let us for example proceed to the late confinement for (15^ 
obtained by requiring that the singularity confine at the level of Xn+i 4 - We obtain the constraint 
an+i 3 an+ibn+iobn+ 9 bn+ 5 K+i = ara+i 2 an+ 7 an+ 2 &ra+i 3 &n+i which, when = a„, reduces to the one 
we obtained in the previous section. Taking again we find the characteristic polynomial 

(fc-I)(Ai2-bl)-bA“-bA®-bA-fc(A9-bA®-bA^-bA^) which factorises into + 

and [{k — 1)(A^ -b I) -b A]. The first of these factors is the one obtained when we implemented late 
confinement on (1551) . More importantly however, the second factor is common to both the standard 
and the late confinement for (ESI) and its role is now clear: it is induced by the special choice of 
gauge, bn = an- Such a gauge-induced factor has no bearing whatsoever on the integrability of the 
mapping and should be disregarded when analysing the characteristic equation. 

Thus, the correct prescription for the implementation as an integrability detector of singularity 
confinement, with full-deautonomisation, can be codified as follows: one should not only perform 
the analysis for the standard confinement but also for the late ones (although, often, the first late 
confinement will suffice). Any common factor in the characteristic polynomials of the standard 
and late confinements indicates the existence of a gauge freedom. As such, it does not play any 
role in the integrability of the mapping and should be disregarded. 

Lest this argument appear arbitrary we proceed now to offer its algebro-geometric justification. 
As before, we interpret the mapping (1591) as a birational map on x = {xn,yn) U {xn,tn) U 
yn) U ^n) (with Sn — 1/ Xnt tn — 1/ Vn^ • 

tpn- P^XP^—^P^XP\ {Xn,yn) ^ {Xn+l,yn+l) = (vn,— -V (61) 

\ / 


This map has two indeterminate points, {sn,tn) = (0,0) and (a:„,y„) = (^0, which can both 
be regularised with a single blow-up. Indeed, the blow-up at the first point 


(-^n, bn ) ^ 



(62) 


immediately yields the mapping yn+i = f^(antn — bn), and the blow-up at the second point 



(63) 


the mapping y„+i = —h 


^ X P^ —^ P^ X P^ 


, both of them well-defined. The inverse map 

Qn bnXr. 


(^n-t-1, yn-t-1) ' t {Xn,yn) — 


yn 


1 ^n+1 


(64) 
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also has two indeterminate points, = (0,0) and {xn+iTyn+i) = (^f^)Oj, 

also be regularised with just one blow-up each. Their respective coordinate charts are: 

{Sn+l,tn+l) <— ( S„+l, U (—, 

\ Sn+1/ V*n+1 / 


which can 


(65) 


and 




bn 


Un+l I Xn+l — 


bn ^n+1 


Un+l 


U 


/ Xn+1 - 

V Vn+l 


1 Vn+l 


( 66 ) 


These blow-ups are depicted in figure [51 together with the resulting exceptional curves. Note that 
'tpn-i ipn share the indeterminate point (sn,tn) = (0,0), for all n, which can therefore be 
resolved by one and the same blow-up (to be performed at each value of n). 




- 1 I- tn = 0 -( I- 

I > 

-= 0 -•- 

Xjj, — 0 S/i — 0 

Figure 6. Diagram showing the base-points in the two blow-ups for the map ipn 
(on the left) and those for the two blow-ups for the inverse map (on the right), 
together with the resulting exceptional curves. 

If we denote the second indeterminate point of '0)0 -1 as 

Pn ■■ (a;„,j/„) = , (67) 

and that of '0n+i as 

Qn+i ■■ (a;„+i,?/„+i) = ^0, , (68) 

then the possible scenarios for regularising the map for all n, can be summarised as: 

'0n+5^ * * * '0n(7n) — (59) 

for some non-negative integer £. 

In particular, the case £ = 0, i.e. ipn(,Pn) = Qn+i, will be seen to correspond to the early 
confinement discussed in section O € = 1 to the standard confinement and the cases .^ > 2 to all 
possible late confinements for this map. In order to analyse these different cases in all generality 
we first define the points 

: (■Sn,yn) = (a,0), T,(^)(/3) : (x„, j/„) = (0,/3), : {xn,t„) = {-f,0), (70) 

=(0,<5), Ti^He):isn,tn) = {0,e), (71) 
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which are related through the action of '0„: 

, (72) 

MTi^HS)) = Ti% (J-^ , . (73) 

We have that Pn = Tn'^ Qn+i = the case of early confinement 

(^ = 0) corresponds to 

Til) ^ fa^\ ^ ^ 

\an-iJ \ a„_i J 

which yields the condition: 

Hn^n—1 Hn+l 


(74) 


, • (75) 

^n—1 ^n+1 

Note that when the requirement (1741) is satisfied, the map i/>„ becomes well-defined after a mere 
three blow-ups. 

When bn = an, he. in the case of the mapping (15511 . this amounts to requiring On = 1 for all n, 
which is indeed the early confinement condition discussed in section [5] On the other hand, if we 
impose &„ = a*, we find the condition 


which can be solved in terms of the characteristic polynomial 

+ \) + \, 


(76) 


(77) 


which is identical to that for the gauge transformation found above. 

For the case t = 1, the condition that needs to be fulfilled for the map to become well-defined is 


pii) ^ ^ I ^ I ^ ^ 


^n—1 


-,(2) f ^n^n—1 




( 5 ) 


n+4 


^n—1 


bn+3bn+2an-l 


^n—1 




(3) f ^n^ri—l 


n+2 


( 1 ) 

n+5 


^n—1 

bn+3bn+2an-l 




( 4 ) 


n+3 




( 2 ) 


n+6 


^n^n—1 

bn+2an-l 
Hn+SHn^n—1 
bn+3bn+2an-l 


— Qn+6 , (78) 


which will only take place after 8 blow-ups since the intermediate points ,■■■ , 7)1+5 ^^1 each 
require a separate blow-up as well. The condition on the parameters that arises from the above 
constraint is of course 

Hn+S^n^n—1 IIn+6 


(79) 

(80) 


bn+3bn+2an-l bn+6 ’ 

which, in case b = , yields 

„fc—1 —1 _ r,^ 

“n+6“™-l-5Una„_i — U„+3Ujj_|_2. 

The characteristic polynomial for this relation is 

(fc - 1)A’^ -b A® - fcA"* - fcA® -f A -b (fc- 1) = (A - 1)2(A -b 1)(A2 -b A -b 1) [{k - 1)(A2 -b 1) -b A] , (81) 

as found above for the standard confinement. The corresponding condition for the mapping (j55p 
is recovered at fc = 1. 

For general i > 2 the closure condition (l69l) takes the form 


(1) / bn-l 
"■ VOn-l 






( 2 ) 


^n-\-5j 


n+5^+1 


^n^n—1 


(^n—1 ^n+5j —2^n+5j—3 


n 


— Qn+l+5^5 


which amounts to the constraint 

^n^n—l 


n 


^n-\-5j 


(^n+d£+l 


On-l bn+5j-2bn+5j-3 bn+5l+l 


(82) 


(83) 
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on the parameters. When we put b = this yields 

e 




^n+Sj 


j = l ®n+5j-2®n+5j-3 


= 1 , 


(84) 


which can be solved in terms of the characteristic polynomial 

e-i 

ge{X) := {k - 1) (A^'^+2 + l) + X + X^{X^ - kX - k)Y^ X^^. (85) 

j=o 

At .^ = 2 we find 32 (A) = (fc — 1)(A^^ +1) + A^^ + A® + A — fc(A® + A® + A^ + A^), which is nothing but 
the characteristic polynomial found above for the first late confinement. Note that regularisation 
of the map ipn under the condition (IMl) will be achieved after exactly 5^ + 3 blow-ups of x . 
The resulting surface, of Dynkin type ^ 4 ^^ is depicted in figure 7. 


D 4 



Figure 7. Representation of the exceptional curves in the surface on which tpn 
(Eli) acts as an automorphism. The labels {x = 0} and {y = 0} refer to the strict 
transforms of the corresponding curves in P^ x P^. 


The Picard group of this surface has rank M+b and we choose Di, D 2 , D 3 , D 4 , D 5 ,Ci,C 2 , ■ ■ ■, 
as its basis. The map ipn induces the following automorphism on the Picard group: 


where the curve is linearly equivalent to 


Di — y D2 —^ 7A3 —^ D4 —^ — y Di 
■ Cl —>■ Csf+i —^ C51+2 ~ 




( 86 ) 

(87) 


— Di -\- D3 -1- iA 4 -|- ^~~^(—C54-1-1 + Csj+a -l- C54-1-4). 
1=0 


( 88 ) 


The sequence (157)) corresponds to a general singularity pattern of length 5£ + 3, as found for the 
general (i.e. early, standard and late) confinements of the mappings (1551) and (1551) . 


The automorphism can be represented by 


corresponds to the permutation of the D curves in 


a 

* 

0 

A 


that 


1, where tr is a 5 x 5 unitary matrix 
(which has period 5) and where the 5£ x 5£ 
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matrix A is given by: 


/O 

1 

0 


A = 


0 

1 0 
0 1 0 
0 1 
0 


0 

1 

1 

0 


-1 

0 

0 1 
1 0 1 

V 0 10/ 

As explained before, it is the characteristic polynomial of this matrix A, 


fiiX) := det(A/ - A) = + (-A^ 


A2 + i)^A5^ 


j=o 


(89) 


(90) 


that will determine the possible integrability or non-integrability of the map Tpn- At £ = 1, we find 
/i(A) = (A + 1)(A^+A + 1)(A — 1)^, which is nothing but the characteristic polynomial for the 
integrable case of (|55|) . Discarding the trivial periodic case i = 0, this is in fact the only integrable 
case for the map '0^ as well. It is easily checked that /^(l) = 1 — which is negative for aW I > 2. 
Hence, for any late confinement, the automorphism will have an eigenvalue that is greater than 
1 and will therefore always have positive algebraic entropy and will thus be non-integrable. Note 
that, for £ = 2, this (largest) eigenvalue is exactly the one we found for the first late confinement 
in section [SJ 

It is clear that for this map ipn, the characteristic polynomial obtained from the automorphism 
on the Picard group does not coincide with the one obtained from the constraints on the param¬ 
eters. However, one can readily check that the characteristic polynomials ge{X) obtained from 
the constraints on the parameters a„ and = On, can be factorized in terms of fe{X) and the 
polynomial (fc — 1) (A^ -I-1) -b A which characterizes the gauge transformation between the mappings 

(15^ and (1551) : 

g,(X) = MX) [ik-l)iX^ + l) + X]. (91) 

Obviously, the automorphism ■0* on the Picard group is completely oblivious to the effect of 
gauge transformations: the curves that make up the Picard group are only defined up to linear 
equivalence and the only effect a gauge transformation will have is that of a dilation of the base 
points of the blow-ups that go into the construction of the surface. On the other hand, because 
of the latter effect, the conditions (15^ or (1551) will feel the influence of the gauge transformation, 
which explains why the polynomials fe{X) are free of these spurious factors whereas the poynomials 
ge{X) are not. This justifies the ‘full-deautonomisation cum late confinement’ approach set forward 
in this paper. However, since any spurious gauge freedom is bound to show up in the form of a 
common factor in the characteristic polynomials for the standard and late confinements (and, if 
this possibility exists, in any early confinement), there is no need to perform a detailed algebro- 
geometric analysis such as that performed here, in order to be able to discard the influence of 
excessive gauge freedom when assesing the integrability or non-integrability of a mapping with the 
criterion we propose. 


7. Late confinement of non-integrable mappings 

In section[5]we have considered the deautonomisation of an integrable mapping and we examined 
the consequences of both early and late confinements, ‘late’ being defined as corresponding to 
singularity patterns that are longer than that for the autonomous case. However, the possibility 
for late confinement does not only exist in the integrable cases: whenever a mapping confines, 
late confinement opportunities always exist. In this section we shall illustrate this fact through a 
few examples, examining along the way the consequences of this late confinement on the algebraic 
entropy. 
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We start with the example used in section [2] to illustrate the full-deautonomisation procedure: 


^n+l ^n—1 



(92) 


Let us for the time being neglect the term b/x and study the confinement associated with the term 
Xjx'^. As we have seen, the standard singularity pattern {0, oo^,0} leads to an equation for a„ of 
the form a„+i —Un-i = 0 and a simple gauge allows us to take an = 1. However, more confinement 
opportunities do exist: the next one corresponds to the pattern {0, oo^, 0, oo^, 0}, which results in an 
equation for a„ of the form 0^+2 —an + a„_2 = 0 ; the one after that will be {0,00'*, 0, oo"', 0, 000}, 
with a„+3 — Un+i + ttn-i — a„_3 = 0 and so on. In order to obtain the algebraic entropy of these 
confining mappings we introduce the term b/x. As we have seen, the standard singularity pattern 
{0, 000} - corresponding to a„ = 1 - leads to an equation for of the form bn+i — 4bn + bn-i = 0 
and the characteristic polynomial Hq (A) = A^ — 4A +1. Now let us suppose that does not satisfy 
this equation and thus that the singularity is not confined after three steps. The next confinement 
opportunity corresponds now to the pattern {0, 00^, 0,00,0, oo"*, 0} and yields an equation for 
of the form bn +3 — 46„+2 + ^n+i — bn + bn-i — 46„_2 + bn -3 = 0. The resulting characteristic 
polynomial is Hq^(A) = A® — 4A^ + A^ — A^ + A^ — 4A + 1 with largest root approximately equal 
to 3.805962. Again, more confinement opportunities do exist, leading to patterns containing an 
arbitrary number of repetitions of the block {0, 00^, 0,00} ending by {0, 00^, 0} and the associated 
equation for 6„. 

As in section [S] we can explicitly obtain the limit of the algebraic entropy for successively 
delayed confinements. The first simple case corresponds to the confinement being delayed because 
the coefficient a„ does not satisfy the lower constraint equations. In this case, starting with the 
characteristic polynomial for the regular confinement Hq (A) = A^ — 4A + I, we find for the m-times 
delayed one, the expression 

= ;^ 2_1 {PoW - 1 ) + 1 . ( 93 ) 


At the limit m —>■ 00 the largest root of is given by that of Po{X) — 1, which is exactly 4. 

The second simple case corresponds to a„ = I, the confinement this time being delayed be¬ 
cause bn does not satisfy the lower constraint equations. Here we find for the fc-times delayed 
characteristic polynomial the expression 

PoW = ^43y(^o (A) - n°(A)) + P^W, (94) 

or an equivalent one noticing that Pq{X) — Pq{X) = A^(APq (A) — 1). Thus at the limit fc —>■ 00 the 
largest root of Po'(A) is given by that of APo(A) — I, approximately equal to 3.8063007. However, 
it is possible to have a situation where the confinement is delayed m times because of the choice 
of an and a further k times because of The corresponding singularity pattern has the following 
structure: the basic block {0, 00is repeated m -I- I times ending by a 0, forming thus what we 
shall call the block-m ; the block-m at the end of which a simple 00 is added is then repeated k 
times ending by a last block-m. The corresponding characteristic polynomial becomes thus 

P^W = a2^+4_i (a^"-+^(AP°(A) - 1)) +p° (A). (95) 


The second example is related to the H-V mapping, already presented in the introduction and 
which was treated within the full-deautonomisation approach in m- In order to illustrate the 
effect of a late confinement we shall work with the form 

Xn+l + Xn-l = Xn -1 J" . (96) 

The standard singularity pattern of (19611 is {0, oo^, oo^, 0} where the meaning of powers in the 
infinities was explained in section [5J The confinement condition on a„ is just a„+3 = an, i.e. a„ is 
periodic with period 3, while for bn we obtain the equation bn+s — 26„+2 — 26„+i -|- = 0. The 

characteristic polynomial of the latter is Pq (A) = A^ — 2A^ — 2A-I-1, the largest root of which is (3-1- 
•\/5)/2 (approximately equal to 2.618034), its logarithm giving the algebraic entropy of the mapping. 
A first late confinement opportunity, when a„ does not satisfy the first confinement condition, is 





SINGULARITY CONFINEMENT AND FULL-DEAUTONOMISATION 


21 


related to the singularity pattern {0, oo^, oo^, 0, oo^, oo^, 0} and leads to two constraints. The one 
on Gn is a„+6 — a„+3 + a„ = 0 which means that a„ can be expressed as a sum of powers of the 
ninth roots of unity and . (Notice that the choice On = 1 is not possible at 

this stage, contrary to the case of the standard H-V mapping). For we obtain the equation 
bn+e — 2 bn +5 — 2 bn +4 + bn +3 — 2 bn +2 — 26„+i + 6„ = 0 and a characteristic polynomial Ti'(A) = 
A® — 2A® — 2A'^ + A® — 2A^ — 2A +1. Its largest root is approximately equal to 2.727069. A numerical 
estimate of the degree growth using Halburd’s Diophantine approximation method gives, after 12 
iterations, a value of 2.727167 converging nicely towards the value of the largest root. 

As in the previous cases the upper bound for the algebraic entropy for an infinitely delayed 
confinement due to a„ can be easily obtained. Starting from the characteristic polynomial Fq (A) = 
A® — 2A^ — 2A +1, for the shortest confinement we obtain for the m-times delayed one the expression 

= ;,3_1 (^o(A) - l) + 1, (97) 

and using the same reasoning as in the previous cases we find that the largest root of Pm{X) for 
m —>■ oo is given by the largest root of To (A) — 1 which is equal to 1 + -x/S- The logarithm of this 
number constitutes the upper bound for the algebraic entropy for the H-V mapping in case of a 
late confinement. However, it is always possible to satisfy the constraint for a„, for instance taking 
a„ = 1, and delay confinement k times by not satisfying the first constraints on bn- In this case 
the singularity pattern comprises a block {0, oo^, oo^, 0, oo, oo} repeated k times ending by a block 
{0, oo^, oo^, 0}. The characteristic polynomial is given by 

^0 (A) = - a - i)a 4 + pO(A). (98) 

For fc = 1 we find that the largest root of Pq^ is approximately equal to 2.67856. When k —>■ oo, 
the largest root is that of the polynomial A^Pq (A) — A — 1 which turns out to be approximately 
equal to 2.678712. 

Our final example is based on the extension of the H-V mapping, introduced in m and analysed 
also in m, 

Xn+l + Xn -1 = Xn -\ -1-(99) 

where we assume that the exponent k is an odd number. As was shown in m a sufficient condition 
for ()98|) to possess the same confined singularity pattern as the standard H-V mapping is that 
a„ = (—1)". When this condition is implemented we find that must satisfy the confinement 
relation bn +3 — k{bn +2 + bn+i) + &« = 0 leading to the characteristic polynomial A® — fc(A^ + A) -I- 1. 
On the other hand, taking a„ = 1 leads to a mapping with unconfined singularities. In m 
Kanki and collaborators have derived the exact value of the largest root of the characteristic 
polynomial for this unconflned case, obtaining the expression (fc + ^Jk(k + 4))/2. It is interesting 
that this value can be easily obtained within the late confinement approach. In fact, an infinitely 
delayed confinement is nothing but a case with unconflned singularities. Following the derivations 
presented above it is easy to convince oneself that the largest root of the characteristic polynomial 
for infinitely delayed confinement is given by the largest root of the polynomial corresponding to the 
shortest confinement minus I. In our case this is the largest root of the equation A^ — A:(A + 1) = 0 
which leads, as expected, to the very same value obtained by Kanki and collaborators. Thus it is 
possible, using singularity confinement, to obtain the algebraic entropy of a mapping which does 
not confine, provided it possesses one special confining case. 

8. Summary and Conclusions 

Singularity confinement is the discrete analogue of the Painleve property. Just as the latter en¬ 
sures that the solution of a differential equation can define a function by eschewing multivaluedness- 
inducing singularities, the former ensures that a mapping is well-defined by resolving all indeter- 
minacies. However, discrete integrability is in some sense more demanding that its continuous 
counterpart since it is intimately related to the growth properties of the solution of a given system. 
Singularity confinement alone might, in principle, be thought to be able to guarantee the slow 
growth that accompanies the integrable character of a rational mapping, since it is based on the 
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simplification of the relevant factors that appear during the iteration of the mapping. This, to¬ 
gether with the observation that all known mappings integrable through spectral methods possess 
confined singularities, was what led two of the present authors, in collaboration with V. Papa- 
georgiou, to advance the hypothesis that the singularity confinement property could be used as an 
integrability detector [I]. 

However, the usefulness of singularity confinement as a discrete integrability criterion was ques¬ 
tioned by the discovery of a counterexample - the (in) famous H-V mapping - that possesses con¬ 
fined singularities despite being non-integrable. This fact cast serious doubt on results obtained 
by singularity confinement, even in the case where the use of the method was perfectly justified, 
for example in the deautonomisation procedure. The aim of the present work, already initiated in 
na, is to show that singularity confinement can indeed, when implemented correctly, be a reliable 
and efficient discrete integrability detector. It goes without saying that in order to do this, one 
must explain the apparent paradox of non-integrable mappings with confined singularities. 

The origin of our present approach can be traced back to the algebro-geometric justification of 
the deautonomisation procedure we presented in m- It was shown there that the constraints on 
the parameters obtained from singularity confinement were equivalent to the linear transformation 
induced, by the mapping, on part of the Picard group of the (family of) rational surfaces on which 
it can be regularised by blowing-up from x P^. This phenomenon does not depend on the 
integrability of the mapping and we were thus led to the conjecture that for any confining mapping 
of the plane, the behaviour of the parameters of the mapping satisfying the confinement constraints, 
is governed by the linear action on the Picard group, obtained after successfully blowing-up the 
mapping. Moreover, as shown in [8], knowledge of the action on the Picard group allows one 
to calculate the algebraic entropy of the mapping rigorously. This observation was only one step 
away from the conjecture that if the deautonomisation of a confining mapping is sufficiently general 
so that its parameters depend on the largest eigenvalue of the linear map on the Picard group, 
the value of this eigenvalue and hence also the algebraic entropy of the mapping can be, so to 
speak, read-off directly from the confinement constraints themselves. The way to implement this 
prescription systematically is through the approach of full-deautonomisation we introduced in m- 
The latter consists in considering a non-autonomous extension of a given mapping, not limited 
to the terms already present, but by considering all the terms that one can add to the mapping 
without modifying its singularity structure. 

All known examples of confining non-integrable mappings, including the H-V mapping, have 
been explained within this new framework for singularity confinement. We were not only able to 
show why one does not expect these mappings to be integrable, but this was done by obtaining 
their algebraic entropy directly from the confinement constraints. This is a great simplification 
since the traditional way of computing the algebraic entropy of a confining mapping is either 
by painstakingly regularising the mapping through successive blow-ups, ultimately obtaining the 
associated action on the Picard group, or by iterating the mapping up to sufficiently high orders so 
that one can establish a recursion relation for the degree growth with some degree of confidence. 

The notions of early and late confinement have been revisited in this paper. Early confinement, 
whenever possible, corresponds to the first possibility to confine and is already present in the 
autonomous case. This early confinement leads either to a trivial, usually linear, mapping or 
to a periodic one. Late confinement, on the other hand, can only exist in the non-autonomous 
case as it corresponds to a singularity pattern longer than that for the autonomous one. Since 
late confinement invariably leads to a non-integrable system, one might of course question its 
usefulness. However, as we have shown, analysing a late confinement along with a normal one can 
help identify any inherent gauge freedom in the mapping, which, when treated in an improper way, 
might otherwise lead to wrong conclusions concerning its integrability. Of course, gauge freedom is 
totally transparent in an algebro-geometric treatment, but it is reassuring to be able to deal with 
this potential difficulty at the level of singularity analysis (after all, not everybody is a superhero 
who is able to find the proper gauge while standing on his head). Since a confining mapping, 
generically, has an infinite number of periodically recurring late confinement opportunities and 
since at each step new confinement constraints will appear, there exist infinitely many possible 
values for the algebraic entropy for a deautonomised mapping, all of which can be simply obtained 
from the singularity confinement approach. 
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Finally, it is our belief that the singularity confinement criterion presented here, based on the full- 
deautonomisation method, can be cast upon a solid, rigorous, basis through an algebro-geometric 
approach and we intend to return to this question in the near future. 
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Appendix 

The eight base-points Qn'' (i = 1 ■ ■ ■ ,8) needed in the regularisation of the indeterminate point 
{sn,yn) = (0,0) for the mapping (j)n (as defined by (15^ 1. together with the coordinate charts 
introduced in each blow-up and the resulting expressions for the mapping: 

i) blow-up at Qn'^: (s„,?/„) = (0,0) ; coordinate chart : 


{Sn, 2/n) ^ ( S„, — ) U ( —, 
/ \ Un 


^ Vn 


Vn+l = ( -Vn + bnVn— + O' 


ii) blow-up at Qn^: = (0,0) ; coordinate chart : 

f Sn \ , f Sn yl\ , , f Sn \ _ 1 f 2 , , 3 

\yn / \yn XUn / ^ 

iii) blow-up at Qn"' '■ (^^,yn'^ = (0,0) ; coordinate chart : 

n \ , f Sn yn\ , , f Sn \ _ 1 f i h 5 Sn Sn 

2 , 2/n 1 ^ ( 2 > ) ^ ( 7,3 > 2^" ) ^ “ ,,4 Sri. ( ^^^^n 3 + “n 3 

n J \yn Sn/ XUn / ^n y3 \ yn yn 

blow-up at Qn'^: (^^,yn'^ = (0,0) ; coordinate chart : 

''n \ , f Sn yt \ , , f Sn \ _ 1 f u 3 Sn , ^ 

3 , 2 /n 1 ^ I 3 , [y4^ynj ^ Vn+l “ 4 ^ 4 + 

'n J \yn / \yn J ^ ^ 

) blow-up at Q^n '^: (^^,yn'^ = J coordinate chart : 


Sn ^ \ ( Sn 1 

4 ^yn J I 4 , 

Vn an J V yt an 


1 


Wn Wn “n 




1 ^ \ , 7 2 

Z/n+l = 3 Sn [ ( ~4 ) 

7/c5 ^ V V 7,4 J 


d) blow-up at Qn^: ,yn'^ = (0,0) ; coordinate chart : 


yn \ Vn an 


1 


,yn] ^ ( — ( ^- 

\yn \yt an 


1 


Vn nn 


U 


1 


1 


2 1 „4 


Vn \yn an 


I i I L 

^ ( TJ ( 71 - “ j +^"2^"F4 

Un y4 \yn \ y-n J Un 


( 100 ) 


( 101 ) 


( 102 ) 


(103) 


(104) 


( 105 ) 
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vii) blow-up at , j/„^ = (0,0) ; coordinate 


1 

vl \y 


77 TJ “ “ ’ 2/" ^ 


1 


vi \y 


yr, 


chart 


1 


U 4 


Vn+l — 


1 


-± \yn \yi 


f '^n ^ \ \ U 

yn^ V^n KVn J ” 2/^ 

hi) blow-up at Q^n'’■ ,2/n) = (-^>o) 5 coordinate 


( 106 ) 


Vlll 


chart : 


1 

Vn \y 


1 \ \ i ^ ^ i I 

71“— +77>2^" ^ 73 71 “— +7T 


VnKyt an) 

\y^ an J 




1 /a„ / i / 6-„ c \ 0„ \ / 6-„ i 

^"+1 = — — 73 74-— +7J +^" 74- — 


fr V2/« \yn \yt 


Vn V2/« Vy'‘ 
b 


Vn 



(107) 


The base-points and coordinate charts needed in the successive blow-ups of the indeterminate 
point {xn+i,tn+i) = (0,0) of ^7^ : 


i) blow-up at Pn)i : (x„+i,t„+i) = (0,0) ; coordinate chart : 

(a;„+i,t„+i) ^ (xn+i, U 

V a;„+i J \ tn+i J 

ii) blow-up at Pnh : (^Xn+i, = (Oj 0) ; coordinate 


^n+15 


^ (a 
Xn+l J \ 


tn+1 


chart : 


re+l in+1 


Xn+U^— U 

Xi+ij \tn+l Xn+l 

iii) blow-up at Pnh : (^Xn+i, = (0, 0) ; coordinate chart : 

f tn+l \ f tn+1 \ f Xn+l tn+1 \ 

\ Xn+l, 2 )^(^"+l’a;3 j ^ i t , 1 ’ a ;2 j 

\ Xn+l/ \ Xn+l/ \ '-n+l X^_^_i / 

iv) blow-up at P^i : {^n+i, = (Oj 0) i coordinate chart : 

( tn+l \ ( tn+1 \ f Xn+l tn+1 \ 

r"+^’ xij ^ xij ’ xij 

v) blow-up at Pn)i : (xn+i, = fo, ; coordinate chart : 


(108) 


(109) 


( 110 ) 


( 111 ) 


^ ^n+1 

^n+15 ~~A 

<+l 




1 ftn+1 1 


Xn+l \^n+l 


u 


Xn+l 


^n + l 

- 


1 ’ ^4 


tn+l 1 

Ojn 


( 112 ) 


'^n+1 


vi) blow-up at : (^Xn+i, “ ^)) = ’ coordinate chart : 
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) blow-up at : {xn+i, ^ ’ coordinate chart : 


1 


^n+l5 2 


^n+1 


.+ 

^n+1 


a + 

-1 1 ^ I Xn+1 


n+1 


^n+1 


^n+l 


^n+1 


u 


^n+1 


-^ X 


n+1 


t'n+l 

^n+1 


- - (114) 


viii) blow-up at : (^Xn+i, “ ^)) = (®’ ’ coordinate chart : 


Xn-\-l 1 


^n+1 


f ^n+1 1 ^ 


\ 




1 


(I'n+l 1 




1 


u 


(^n+1 


n+1 

1 


^n+1 

^n+l 


^n+l 


hi \*^n + l / “n 



(115) 
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